On Koopman-von Neumann Waves 



D. MauroQ 

Dipartimento di Fisica Teorica, Universita di Trieste, 
Strada Costiera 11, P.O.Box 586, Trieste, Italy 
and INFN, Sezione di Trieste. 



In this paper we study the classical Hilbert space introduced by Koopman and von 
Neumann in their operatorial formulation of classical mechanics. In particular we 
show that the states of this Hilbert space do not spread, differently than what happens 
in quantum mechanics. The role of the phases associated to these classical "wave 
functions" is analyzed in details. In this framework we also perform the analog of the 
two-slit interference experiment and compare it with the quantum case. 



e-mail: mauro@ts.infn.it 



1 Introduction 



In their standard formulation classical and quantum mechanics are written in two 
completely different mathematical languages: for example in classical mechanics ob- 
servables are functions of a 2n-dimensional phase space, while in quantum mechanics 
they are self-adjoint operators acting on an Hilbert space. In the literature there are 
a lot of attempts to reformulate classical and quantum mechanics in similar forms, 
|l|]-[@]. In this paper we shall concentrate on the work of Koopman and von Neumann 
(KvN) who proposed, in 1931-32, an operatorial formulation of classical mechanics, 
^-0]. The starting point of their work is the possibility of defining an Hilbert space of 
complex and square integrable classical "wave" functions V'(v^) such that p{ip) = \ip{ip)\'^ 
can be interpreted as a probability density of finding a particle at the point ip = (g, p) 
of the phase space. This p has to evolve in time according to the well-known Liouville 
equation: 

d 

^-g^p{(l,p) ='Hp{q,p) (1.1) 

where 7i is the Liouville operator Ti. = —idpHdq + idqHdp and H is the Hamiltonian 
of the standard phase space. In order to obtain ( |1 . 1|) Koopman and von Neumann 
postulated the same evolution for ifj: 

d 

i-Q-^i^{<l,p) =T<i){q,p) (1.2) 
Since Ti contains only first order derivatives it is easy to check that from eq. ([1.2D one 



can derive eq. (|1 . 1| ) . This is something which does not happen in quantum mechanics 
where the analogue of Ti is the Schrodinger operator which contains second order 
derivatives. 

Now in quantum mechanics the complex character of the wave function is of funda- 
mental importance for a lot of reasons: while the modulus of the wave function gives 
the probability density p, the phase of ip brings in also some physical information. In 
fact it is related to the mean value of the momentum operator p and it gives origin to 
the appearance of interference effects in two-slit experiments. To answer the question 
if phases of classical "wave" functions ^(g,p) play an analogous role also in classical 
mechanics is the main goal of this paper which is organized as follows. 

In section 2 we will briefly review the Koopman-von Neumann formalism and the 
associated functional formulation 0. In section 3 we will present a very simple but 
pedagogical example that shows how, in the operatorial approach to classical mechan- 
ics, the lack of an uncertainty relation between q and p and the different form of 7i 
with respect to the Schrodinger operator H prevents the spreading of the wave func- 
tions. We shall also show that the phase and the modulus totally decouple in the 
equation of motion of ip and that the phases do not influence the expectation values of 
the observables of classical mechanics. In section 4, after spending some words about 
the abstract Hilbert space of classical mechanics, we shall underline how a lot of the 
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previous considerations are a consequence of the particular representation we have cho- 
sen. It is the representation where the classical "wave" functions ip{q,p) are given by 
functions of q and p. If we change representation the situation changes drastically. In 
particular if we use a representation where p in ijj{q, p) is replaced by a different variable 
Ap, then phases begin to play a crucial role since they bring in physical information 
(e.g. the mean value of p). In this representation Ap) and p{q, Ap) do not evolve in 
the same way, as ip{q,p) and p{q,p) used to do, and there is, in a certain sense, a sort 
of spreading of the wave function. Since in this representation phases seem to play a 
physical role we can say that it is necessary to consider, for mathematical consistency, 
an Hilbert space made up with complex and not with real wave functions. 

In quantum mechanics the complexity of the wave functions is one of the most 
important reasons for the interference effects. Therefore a question that arises quite 
naturally is: what happens in this operatorial approach to classical mechanics if we 
consider complex wave functions? In section 5 we shall propose the classical analog 
of the two-slit experiment and we will perform in details all the relevant calculations. 
We will show that the different form of the evolution operators in the quantum and in 
the classical case leads to totally different results: the known interference phenomenon 
at the quantum level and its non-appearance at the classical one. As we said above 
we tentatively link this difference to the different form of the evolution operators but 
more work has to be done to get a deeper understanding of the phenomenon. 



2 Operatorial Approach to Classical Mechanics 

In quantum mechanics, starting from the Schrodinger equation for the wave function 

ih^^^m^,t) =, in^^ = -^^i.{x,t) + V{xmx,t) (2.1) 

we have that the probability density p{x,t) — |'0(x, satisfies a continuity equation 
of the form: 

^ = -dtvj (2.2) 
where we have indicated with j the probability density current: 

j^-^(^r^^-i^^r) (2.3) 

Now if we write the wave function as ip = y/pexp[iS/h] we discover immediately that 

— * 

the phase S enters explicitly into the expression of the current probability density j: 

J=^ (2.4) 

m 
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As a consequence the equation of evolution of p ( p.2| ), couples the phase and the 
modulus square p of the wave functions, 0. We can also notice from eq. (p.2|) that in 
quantum mechanics the probability density p does not evolve in time with the same 
Schrodinger Hamiltonian H which gives the evolution of the wave function if). 

The situation is completely different in the Hilbert space of classical mechanics. In 
fact, following Koopman and von Neumann, we postulate that the wave functions 
ip{ip^t) = ip{q-,P-,t) evolve in time with the Liouvillian operator: 

n = -id,m,^ + id,^Hd,^ (2.5) 

according to the following equation: 

^^^i^ = n^ ^^i^ = {-d,^Hd,^ + d,m,:)^ (2.6) 



We can think of ( p.6| ) as the analogue of the quantum Schrodinger equation, i.e. as 
the fundamental equation governing the evolution of the vectors in the Hilbert space 
of classical mechanics. These vectors are the complex wave functions on the phase 
space obeying the normalizability condition / dqdpip*{q,p)ip{q,p) = 1. If we take the 
complex conjugate of ( p.6| ) we obtain: 

^^^* = {-d,^Hd,^ + d,,^Hd,^)r (2.7) 



i.e. ip and ip* satisfy the same equation. Now if we multiply (2^) by ip* and (2.7) 



by ip cind we sum the two equations we re-obtain eq. ( |1.1| ), i.e. the evolution of 
p{q,p) = ip*{q,p)il'{q,p) with the Liouvillian operator: 

-p=(-d,^Hd,^ + d,^Hd,^)p ^ ^^P = ^P (2-8) 



So we have derived the standard Liouville equation ( p.8|) as a consequence of having 
postulated eq. (^.6]). Moreover we notice that ( ^.81 ) does not couple the modulus square 
of ip, i.e. p, with the phase of ip differently from what happens in quantum mechanics, 
eq. (|23) and (0). 



If we define the scalar product between two wave functions and r as {iP\t) = 
J d^pip*{ip)T{ip) it is easy to show that {i/jIHt) = {Hiplr), i.e. the Liouvillian 7i is a 
self-adjoint operator; consequently the norm of the state {ip\ip) = J dip'ip*{ip)ip{ip) is 
conserved during the evolution and we can consistently interpret p{ip) = ip*{ip)ip{{p) as 
a probability density in the phase space. 

Before going on we want to show how the Liouvillian 7i arises in a natural way in 
the functional approach to classical mechanics described in ref. which represents 
the path integral counterpart of the KvN operatorial formulation. Let us ask ourselves: 
which is the probability of finding a particle at a point y?" = (g\ . . . , g"; pi, . . . ,p„) in 
phase space at time t if it was at (pi at the initial time tj? This probability is one if 



3 



(f^ and are connected with a classical path 0"^, i.e. a path that solves the classical 
equations of motion, and it is zero in all the other cases. So we can write: 

p(^^#^t,,) = %'^-C/(^;^.)) (2.9) 

where the RHS is a functional delta that forces us to stay on the classical path 
associated with the initial condition ipi. The functional delta can be rewritten as a 
delta on the Hamiltonian equations of motion ip"' = uj"'^dbH[(p) via the introduction of 
a suitable functional determinant: 

5(<^'^ - C,(t; ^i)) = 5(0" - uj-%H)det{dt6^ - u'^'d.d.H) (2.10) 

It can be shown, see eq. (3.51) of ||^, that the determinant in the previous equation 
is formally independent of the fields and so it can be put equal to one if we are 
not interested in the study of nearby trajectories. Now if we exponentiate the Dirac 
delta in ( p.lOp via the introduction of auxiliary variables \a we get the following path 
integral: 

P((^,t|(^i,ti) = J V^di^" - uj^'^dbH) = JVipVXe'I^'^' (2.11) 
where the Lagrangian C is: 

C = Xa^''-n (2.12) 

with 

n = XaUJ^'^dhH (2.13) 

From the kinetic part of the Lagrangian we can deduce the form of the commutators 
of the associated operatorial theory: 

[rX] = ^S^ r,\] = [p\ApJ = z5} (2.14) 

The previous commutators^ can be realized considering as multiplicative operators 
and Xb as derivative ones: 



From now on we will indicate this representation as the Schrddinger representation of 
classical mechanics. Using ( 2.15 ) Ti becomes the following operator: 

n = -tLu'^'dbHda = -idp^Hd,^ + id.^Hdp^ (2.16) 



that is exactly the Liouvillian 7i of eq. (|2.5|) . This confirms that the path integral 
( p.ll| ) is the correct functional counterpart of the KvN operatorial theory. 



^AU the commutators different from (2.14) are identically zero. In particular, differently from the 
quantum case, we have that [q,p\ — which implies that we can determine with an arbitrary precision 
the position and the momentum of a classical particle, like it happens in the standard phase space 
approach to classical mechanics. 
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In rederiving ( 2.11 ) we have started from the transition probability P{(p,t\(pi,ti) of 
going from ifi to in a time interval t — ti. So we can say that the path integral 
( p. lip gives a kernel of evolution for the classical probability densities p{ip,t), in the 
sense that if we know the probability density at the initial time ti we can derive the 
probability density p at any later time t via the standard relation: 

p{ip,t) = dipiP{<^,t\<^i,ti)p{<^i,ti) (2.17) 



Let us remember that KvN postulated for ijj the same equation of motion (|2.6| ) as for 
p. As a consequence their evolution can be represented as: 

^lj{(p,t) = J d(piK{(p,t\(pi,ti)tjj{(pi,ti) (2.18) 

where the kernel of evolution K{(f,t\(fi,ti) has the same expression as the kernel of 
evolution P{ip,t\ipi,ti) for the densities p. The reader may be puzzled that the same 
kernel propagates both ip and {ipl"^- We will see that there is no contradiction in this 
by working out in details the case of a free particle in appendix A. 



3 Spreading and Phases of the Wave Functions 

One of the most characteristic effects of quantum mechanics is the spreading of the 
wave functions during their time evolution. Let us consider a quantum system of a free 

- 2 

particle in one dimension with Hamiltonian H = —Ti —— and initially described by a 
gaussian wave function: 

1 f i \ 

%l){x) = —^=exp{-— + -pix] (3.1) 



2a2 h 



It is easy to check that at the beginning the mean value of the position x is equal to 

zero and the uncertainty in the measurement of the position is: (Ax)^ = — . At time 
t the wave function will be: 



'?/'(x, t) = N ■ exp 



. X 1 a; + T \ Pi^ ^ 

2[ma'^ + iTit) \ mj n\ 2m. 



m ( p%t\^ i ( p^ 



(3.2) 



We can note how the coefficient Pi, which entered only into the phase of the wave 
function at time t = 0, managed to enter also the modulus of the wave function at 
time t > 0, see eq. (|3.2|). As a consequence the expectation value of the position x at 
time t depends explicitly on the original phase p^; in fact we have: 



x{t) = dx X ■ \ilj{x,t)\ =pit/m (3.3) 
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So we see that the information about the mean value of x is carried by terms appearing 
in the original phase of ip. For the mean square deviation of x we get: 

(A^=(x-x(t)P = ^(l + ^) (3.4) 
So, at t — > oo the wave function will be totally delocalized: 



lim {Ax{t)y = +00, Va > (3.5) 

t — >oo 

This effect is present also if we prepare an initial state very sharply peaked around the 
origin, in fact: 

lim A(x(t))2 = limf^ + 7-^) = +00, Vt > (3.6) 

Note that the previous limit is +00 because of the presence of the parameter 0? in 
the denominator. The relation ( ^.6] ) is not surprising: if a — at the beginning then 
we have a state perfectly localized in space, i.e. (Ax)^ and, from the Heisenberg 
uncertainty relations, we can deduce that (Ap)^ — * +00. So in this case the initial 
momentum is completely undetermined and, consequently, even after an infinitesimal 
time interval, also the position of the particle becomes completely undetermined. These 
are the well-known quantum mechanical effects. 

What happens in the operatorial version of classical mechanics? As we have seen 
in the previous section, the evolution in time of the wave functions is generated by 
the Liouvillian itself Ti = —idpHdq + idqHdp which, in the particular case of a one- 
dimensional free particle, has the following simplified form: 

H = (3.7) 

The free Liouvillian is essentially the product of two commutative operators: an opera- 

d 

tor of multiplication p and an operator of derivation —Itt- So if we want to diagonalize 

oq 

d 

the Ti of eq. ( p.7| ) we have to diagonahze simultaneously both p and The eigen- 

states of p associated to an arbitrary real eigenvalue po are the Dirac deltas 5{p — po); 

d 1 
the eigenstates of —i-^ are instead the plane waves exp[iXqq] and their correspon- 



dq 

dent eigenvalue^ are \q. So the eigenstates of the free Liouvillian ( |3.7|) are just the 

d 

product of the eigenstates of p and — ^7:— : 

oq 



VZtt 



d d 

^ We call Xq the eigenvalues of —i-^ since — i-r- is just a representation of the abstract Hilbert 

aq oq 

space operator A^: see (2.15) and the next section. 
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and the associated eigenvalues are the product of the eigenvalues: S = —2^. Now, 

m 

suppose to take as initial wave function the following double gaussian in q and p: 

where a and h are related to our initial uncertainty in the knowledge of q and p: 
(Ag)2 = (Ap)^ = 6^/2. Note that, since in classical mechanics q and p commute, 

there isn't any uncertainty relation. As a consequence a and h in (^.91) are two com- 



pletely independent paramaters and the product (Ag)^ • (Ap)^ can assume arbitrary 
small values. 

Now we can write the initial wave function ( p.9|) as a superposition of the eigenstates 
(B.8|) of the Hamiltonian Ti as: 



i/j{q,p,t = 0) = y dXgdpoc{Xg,po)Tx^p,-,{q,p) (3.10) 
where the coefficients c(Aq,po) are: 

e,p( - (3.11) 



i6 'V 2 262 
The wave function at t is given by: 



'0(9, P^t) = I dXqdpo c(Ag, po) exp[-i£t] Tx^p^ (g, p) 
1 



y/irab 



exp 



2a2 V m J 262 



(3.12) 



Note that this ip{t) is related to ip{0) by the following equation: 

/ vt \ / OH OH \ 
i^{q,P,t)=^[q-^,p,Oj =^[q- -Q^t^P+^t^^) (3-13) 

The previous relation could be inferred also from the path integral ( |2.11|) . In fact in 
the case of a free particle, see appendix A, the kernel of propagation is correctly given 
by: 

K{ip, t\ip,, 0) = 5 fg - - —)6ip - p,) (3.14) 
from which we obtain immediately: 

ilj{(p,t) = J d(piK{(p,t\ipi,0)ilj{ipi,0) = ilj(^q -^,p,0^ (3.15) 
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If we identify the modulus square of the wave function \ip{q,p,t)\'^ = p{q,p,t) with the 
probabihty density of finding a particle in a certain point of the phase space we have 
that ( p.l3| ) implies: 

piQ^P^t) = pyQ- -Q^t^P+ -Q^t^^J (3-16) 

This is an equation in perfect agreement with the Liouville theorem — p = 0. 

at 

Let us now go back to eq. ( p.l2| ) and calculate the mean values of the dynamical 
variables at a generic time t. They are: 

q= I dqdpq-\il){q,p,t)\^ = p = ( dqdpp ■ \il}{q,p,t)\^ = pi (3.17) 

J m J 

Note that, differently from quantum mechanics, the information on the mean values of 
p is given by coefficients which appear in the modulus of the wave function. The mean 
square deviations are: 

2 l2 +2 u2 



¥ t 



(^?W)' = y + Y;;^. (ApW)' = Y (3-18) 



Therefore also in classical mechanics if 6 7^ we have lim (Ag(t))^ = +00 and the 

t — ^00 

wave function is totally delocalized. This is not strange if we consider that we are 
giving a statistical description of a set of particles with a momentum distributed in a 
gaussian way around Pi. This means that we can have particles with momenta both 
greater and smaller than pi. These particles cause a dispersion in the wave function 
and, consequently, in the distribution of probability around the mean value of q. 

If we instead consider the mechanics of a single particle we can measure exactly its 
position and momentum at the initial time. In this case the terms that parametrize 
the gaussians (|3l9| ) go to zero (a — 0, 6 — > 0) and the initial wave function is a good 
approximation for the double Dirac delta S{q)6{p — pi). At the generic time t the state 

will be described by another couple of Dirac deltas: s( q~'^]6(p—pi). In this limiting 

V m J 

case the variances are identically zero because, differently from the quantum case ( |3.6| ), 
in (|3.18|) the parameters and 6^ do not appear in the denominator: 



lim (Ag(t))2= lim +-— =0, lim (Ap(t))2 = lim - = (3.19) 

a,b^O a,b^O \ Z Z Tfl J a,fe^O 6-+0 Z 

and the particle remains perfectly localized in the phase space at every instant of time 
t. 

We feel that even this very simple and pedagogical example can be used to underline 
some very important differences between the quantum and the classical operatorial 
approaches which are: 

1) in classical mechanics we can prepare an initial wave function that approximates a 
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double Dirac delta in q and p since q and p are commuting operators and so there is 
no uncertainty relation between them; 

2) the classical dynamics given by Ti is such that, if we know with absolute precision 
the position and the momentum at t = 0, they remain perfectly determined at every 
instant of time t and there is not any spreading, see eq. ( |3.19| ); 

3) the knowledge about the average momentum of the classical particle is brought by 
terms appearing in the modulus of the wave function. 

For a classical free particle it is easy to show that if we add a phase factor to an 
arbitrary initial wave function this phase factor does not pass into the real part during 
the evolution differently from what happens in quantum mechanics, see eq. ( p.2|) . This 
can be proved as follows: every initial classical wave function 

ilj{q,p) = F{q,p)exp[iG{q,p)] (3.20) 

can always be written as a superposition of the eigenstates of the free Liouvillian ( p.7|) 
in the following way: 

^{q,p) = J dXgdpoc{\g,po)rx^p^{q,p) (3.21) 

where the eigenstates Tx^,po{q,p) are given by eq. ( |3.8| ). So the coefficients c{\q,po) are 
basically the Fourier transform, in the q variable only, of ip{q-,Po)'- 

c(A„ po) = / dq e-'^''^F(g, po)e^''^^'^«) (3.22) 



/27r 

At this point we can find the free evolution of the wave function in the usual way: 



'4'iq,P,'t) = / d\gdpoc{\g,po)exp[-i£t]Tx^p^,iq,p) 



dXq c{Xq,p)exp 



iXq{q-—) 



r [q ,p ]exp 

m 



iG[ q - —,p 
m 



(3.23) 



and we obtain a result that is again in perfect agreement with the kernel of evolution 
( p.l4| ). Therefore, in the case of a free particle, we have that for every initial wave 



function ip{q,p,t = 0) = F{q,p)exp[iG{q,p)] the probability density \'ip{q,p,t)\'^ does 
not depend on the phase G{q,p) not only at the beginning, but also at any later 



time; in fact from (|3.23|) we have that \ip{q,p,t)\'^ = F^iq — —,p]. This has some 



m 

consequences also on the expectation values of the observables. If we assume that 
observables in classical mechanics are only the functions 0{(p) or, in operatorial terms, 
only the operators 0{(p) then it is easy to check that their expectation values do not 
depend on the phase G of the wave function ( |3.2CI| ) : 



(O) = / dipF*{ifi)exp[-iG{ip)]0{ip)F{ip)exp[iG{ip)] = / dipF*{ifi)0{ip)F{ifi) (3.24) 
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This is true also at later times t because of the form (|3.23|) of the wave function at 
t > 0. 

This independence from the phase G would not happen if the observables were 
dependent also on A: 0(<^, A), because A is a derivative operator and it would give 

{O) = |rf^F*(^)ea;pHG(^)]o(^,-z^)F(^)ea;p[zG(^)] = 

= Jd^F*{^)Q{^,F,F',G,G')F{^) (3.25) 

So {O) would be a function of F{ip), its derivative F' with respect to ip, but also of 
the phase G{!f) and its derivative G'. These considerations can be extended from the 
free particle case to a physical system characterized by a generic Hamiltonian 7i. In 
fact the solution of the equation: 

i^^^p{q,p,t)=n^{q,p,t) (3.26) 

is given by ||^: 

il){q,p,t) = tlj{q{q,p,t),p{q,p,t)) (3.27) 
where q and p are the solutions of the equations: 

-/ .^ dH{q,p) ^ dH{q,p) 

q^{q,p,t) = Pj{q,p,t)= ^_ (3.28) 

with the initial conditions qj{q,p,0) = Pj{q,P,0) = p^. So, according to eq. 

( ^.27| ), the evolution of a classical system via the Liouvillian does not modify, in the 
Schrodinger representation, the functional form of ip provided we write it using the 
g, p. This has, as an immediate consequence, that if we take a wave function without 
any phase at the initial time, phases cannot be generated during the evolution: 

H : ipiwithout phases t = 0) — ^ ipiwithout phases t) (3.29) 

In quantum mechanics instead, even if we start from a wave function that does not 
present phases, these ones will be created in general at later times via the operator H: 

H : ip{without phases t = 0) — > ip{with phases t) (3.30) 



This can be seen for example from (|3.1| ) and (|3.2| ). Even if we start with no phase 
Pi = 0, at time t we get that the wave function (|3.2|) becomes: 



ip = N ■ exp 



(3.31) 



2(ma^ + iht) 

and it has a phase because of the term iht in the denominator. All this is a consequence 
of the fact that the phase and the modulus of a wave function interact in quantum 
mechanics as one can see from standard text books, p. In fact writing 



il){x) = A{x)exp 



(3.32) 
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and equating the real and imaginary part of the Schrodinger equation (|2.1|) we obtain 
the following two equations for A{x) and S{x): 

+ 7^^] +V 



dt 2m\dxJ 2mA dx'^ 33") 

^dA dAd^ 4^_o ' 
dt dx dx 2 dx^ 



From (|3.33| ) we see that S and A are coupled by their equations of motion. 



In classical mechanics if we start from 

ip{q, p) = F{q, p)exp[iG{q, p)] (3.34) 

we can insert it in eq. ( p^) : i— = Hip and, equating the real and imaginary part, we 
get: 

dF dC 

= HF, = KG (3.35) 

dt ' dt ^ ^ 

So we see that in classical mechanics the phase and the modulus decouple from each 

other. I owned this analysis to E. Gozzi who has once summarized all this with 

the sentence: " What is quantum mechanics? Quantum mechanics is the theory of the 

interaction of a phase with a modulus" . 

As we have just proved, in classical mechanics there is at least one representation in 

which this interaction is completely lost and the evolution of the modulus is completely 

decoupled from the evolution of the phase. One may then think that it is useless to deal 

with complex wave functions if their phases do not bring in any physical information. 

This is true only if we decide to work with the Schrodinger representation. If, instead, 

we change representation, as we will do in the next section, and use the one where p is 

realized as the derivative with respect to Xp, we shall then show that the mean value 

of p is related to the phase of the wave functions. So, if we want to be as general as 

possible and not just stick to the Schrodinger representation, we have to assume that 

the classical wave functions are complex objects. 



4 Abstract Hilbert space and (g, A^) representation 

In the previous section we have restricted ourselves to the Schrodinger representation 
in which both q and p are realized as multiplicative operators, and we have worked out 
everything in this frame. What we want to do now is to construct the Hilbert space 
of classical mechanics from an abstract point of view independently of any particular 
representation. We can start observing that q and p can be considered as a complete 
set of commuting operators whose real eigenvalues can vary with continuity from — oo 
to +oo: 

q\q,p) = q\q,p); p\q,p) = p\q,p) (4.1) 
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The eigenstates \q.,p) form an orthonormal and complete set which can be used as 
a basis for our Hilbert space. The orthonormahty and completeness relations are 
respectively given by: 



dqdp\q,p){q,p\ 



(4.2) 



The connection between the abstract vectors {ip) and the wave functions ip{q-,p) is given 
by the relation {q,p\ip) = ip{q,p). In this basis the operators q and p are diagonal: 



{q',p'\q\q",P") = Q'5{q'-q")5{p'-p"y. 
{q\p\p\q\p") = p6{q' - q")S{p - p") 



, .1 , .df.d 
while the operators —t^rl "^"^ 

oq V op 



defined by the relations: 



{q',p' 



9_/ ^\ _ _-_^f^_^\>f ^, 
dq\ dp) dq'\ dp' J 



dq' \ dp' 



(4.3) 



(4.4) 



are self-adjoint. From ( [4.2D -( ^74D it is easy to check that: 



{q\p' 



dq. 



^) = {q',p'\i\ 



{q',p' 



. d 

dp 



i)) = {q',p'\i\ 



(4.5) 



while all other commutators are zero. Because of the completeness of {q',p'\ and 
the arbitrariness of we have that ( |4.5|) can be turned into the purely operatorial 

■ i and so we can identify, from eq. ( p.l4|) , 



relations: 



d 



q, 



dq 



i and 



d 



P, 



dp 



Xg = —i— and A„ 
dq 



<9 , ^ d 

-2—. Now it is easy to show that the Liouville equation ( 
dp 

nothing else than a particular representation of the abstract Liouville equation: 

dt 



IS 



(4.6) 



obtained using as basis the eigenfunctions of q and p. If we consider as Hamiltonian 

in the standard phase space the following one: H = h V{q), then the Liouville 

equation (|4.6|) becomes: 



P 



m 



XpdqV{q) 



Projecting the previous equation onto the basis (^,^1 we easily obtain: 



d 

i-g^{q,p\i^,t) = {q,p 



A, 



P 



m 



i^,t) - {q,p\XpdgV{q)\^,t) 



p d d 
~i—Tr{q,p\i^,t)+^dgV{q) — {q,p\ip,t) 
m dq dp 



(4.7) 



(4.8) 
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that is equivalent to the usual Liouville equation: 



dt 



ip{q,p,t) 



m oq op 



ip{q,p,t) 



(4.9) 



The \q, p) basis is not the only one for the Hilbert space of classical mechanics. A very 
important representation [|lOl is the one in which we consider as basis the simultaneous 



eigenstates of q and Ap which, according to ( |2.14| ), are commuting operators: 

g|g, Ap) = g|g, Ap); Ap|g, Ap) = Ap|g, Ap) (4.10) 



Sandwiching the second relation in ( ^.10| ) with the bra {q\p'\ we obtain 

d 

- i-Q^{q',p'\q, Ap) = Xp{q',p'\q, Ap) 



The solution of this differential equation is: 



{q',p'\q, Ap) 



1 



/27r 



S{q - q')e 



(4.11) 



(4.12) 



The states \q,Xp) form a complete set of orthonormal eigenstates, i.e. an alternative 
basis for the vectors of our classical Hilbert space. In this basis we have: 



{q, K 



which, via ( |4.12| ), gives: 



j'dpiq, \p\q',p) {q',p\ip) 



dpe ^P^''^{q,p) 



(4.13) 



(4.14) 



i.e. the wave functions in the new basis are related to the ones in the Schrodinger 
representation by means of a Fourier transform^. In this new representation we have 
for the p operator: 



(g, Xp\p\ip) 



dq'dp'iq, \p\p\q' ,p') {q' ,p'^ 



1 d 
i- 



2tt dXr 



27r 
d 



dp'p'e-'P'^^^lj{q,p') 



dp'e-^^'^^{q,p'\^lj) = i — {qAp\ 



(4.15) 



while for Ap we have: 



(g, Ap|Ap| 



(4.16) 

Summarizing (|4.10| )- (|4.16| ), we can say that in this representation we have to consider 

d 

p as a derivative operator: p = "j^r— and Ap as a multiplicative one. This is simply 



*We indicate the wave functions in the new basis with the same symbol ip for notational simplicity. 
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a different realization of the usual commutation relation: [p, Ap] = i. Using {q,Xp\ to 
sandwich eq. ([4.7|) we get that the Louville equation becomes: 



1 d d 
m dq d\ 



ij{q, Xp) - XpdgV{q)ip{q, Xp) 



(4.17) 



We shall now show that a lot of the results of the previous section were in a certain 
sense representation-dependent. In fact in the new representation, since the momentum 
p has become an operator of derivation, we have that the information about its mean 
value is brought in by the phase of the wave function similarly to what happens in 
quantum mechanics. For example the double gaussian state of eq. ( |3.9| ) becomes the 
following one in the new basis: 



ip{q, Xp,t = 0) 




exp\ 



2a2 



ipiXp 



(4.18) 



We obtain it by just applying formula (|4.14|) and, after the Fourier transform, the wave 
function which was real in the Schrodinger representation becomes complex. The mean 
values of q and p are obviously the same as before: 



\q\ 



0, 



p 



\p\tlj) = {ip 



d 



dXr 



(4.19) 



but now we see that elements appearing in the phases of the wave functions, like pi in 
( [4.18| ), begin to play an important role since they are linked with the mean values of 
physical observables like p. 

Let us now make the evolution of ( |4.18| ) under the Liouvillian of a free particle. This 
Liouvillian in the new representation is given by: 



n 



m dqdXp 

Its eigenstates associated with the eigenvalues pXq/m are: 



(4.20) 



^A,,p(g, Xp) = —exp[iXgq - iXpp] 



(4.21) 



Expanding then the wave function ( |4.18| ) in terms of the Tx^^p above and making the 
evolution of the system, we obtain at time t: 



ip{q, Xp, t) = N ■ exp 



262 



1 {XpUia b + iqtb + ipima 
2 



2\2 



(4.22) 



2a2 262 2 0252(^2^2 + ^252) 

From the previous formula we see how the factor pi which at time t = entered only 
the phase factor, see (|4.18|) , has passed at time t also into the real part of the wave 
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function, exactly as in the quantum case we studied before. The expectation values 
and the variances of q and p are still given by eqs. (p.l7|) - (|3.18|) : 



q = f dqdXp q ■ {^1' = — , P = P^ (4-23) 
J m 

and: 



a2 62 ^2 ^2 



(Ag)2 = - + - — , (Ap)2 = - (4.24) 



This is so since they are observable quantities and, consequently, they are independent 
of the representation we are using. In the new representation making a wave function of 
the form (^4.18|) well-localized both in q and in Xp means sending a ^ 0,b oo. In this 



limiting case we have that (Ag)^ — > cxd at every instant of time t > and so there is a 
sort of spreading of the wave function. This is not surprising. In fact if the initial wave 
function is very peaked around q = Xp = then we know precisely the values of q and 
Xp, instead of the values of q and p. From the commutator [p, Xp] = i we can derive the 
following uncertainty relation Ap- AXp > 1/2, where Ap and AXp are the square roots 
of the mean square deviations. So, if we determine with absolute precision Ap, as we 
do with the limit 6 — > oo, the momentum p is completely undetermined. Consequently 
also the position q at every instant t > is completely undetermined, because it follows 
the classical equations of motion q = p/m. This fact has, as an immediate consequence, 
the spreading of q and the complete delocalization of the wave function at every instant 
following the initial one.0 

Another aspect that we can study is the continuity equation. We have seen in the 
second section that the continuity equation in the Schrodinger representation is nothing 
else than the usual Liouville equation for the probability density p. What happens in 
the other representation we have studied in eqs. ([4.1CI|) - (|4.16|) ? According to what we 
have already seen, the Liouville equation for ip is: 

igli^il, = -^^^(^' ^f) - KV'miQ, K) (4.25) 
while the one for the 4'*{q, Xp) is the complex conjugate: 

- t^^riq, K) = -7r7^^*(^' - ^p^'(l)'^*(l^ ^p) (4-26) 
ot m oq oXp 



From (|4.25| ) and (|4.26| ) we can obtain the equation for p(g, Ap) = '?/'*(g, Xp)ip{q, Xp). It 
is of the form: 

^p(g, Ap) + J = (4.27) 

where: 

''The usual mechanics of the single particle can be reproduced also in this representation but we 
have to take the limit a ^ 0, 6 ^ 0, i.e. we have to use for Ap a plane wave of the type exp{—ipiXp). 
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So we notice that in this case p(g, Xp) evolves with an equation that is completely 
different from the Liouville equation and there is no manner to write J in terms only 
of p. Moreover, if we write ip{q,Xp) = y/pexp[iS{q, Xp)], the phase S{q,Xp) will enter 
explicitly J and the equation of p, i.e. we have a situation very similar to that of 
quantum mechanics where phases and modulus are coupled in the equations of motion. 

Another aspect that the (g, Xp) representation of classical mechanics and the standard 
quantum one have in common is that, even if we prepare a real wave function of q and 
Xp at the initial time t = 0, phases will be created in general by 7i during the evolution. 
This can be seen by means of our usual example. In fact, if we put pi = 0, we have 
from eqs. ( ^A8D and (^^) that: 




exp 



2a? 



N ■ exp 



q^ 1 {XpTna'^b^ + iqtb 



2\2 



2a? 2aW{m?a^ + m 



(4.29) 



that is: 



Ti : ip{q,Xp,t = Q without phases) 



4'{q,Xp,t with phases) 



(4.30) 



Since the (g, Xp) representation has all these features in common with quantum me- 
chanics it is not a case that this representation turns out to be [|10] the one where the 
transition to quantum mechanics is best understood. 



5 Two-slit experiment 

Having formulated classical mechanics in the same mathematical language of quantum 
mechanics, we think it may be useful to analyze the two-slit experiment in the classical 
KvN formalism and compare it with its quantum analogue. This kind of experiment 
is central in quantum mechanics and its mystery is best summarized in these words of 
Feynman [^: "The question is, how does [the two-slit experiment] really work? What 



machinery is actually producing this thing? Nobody knows any machinery. The math- 
ematics can be made more precise; you can mention that they are complex numbers, 
and a couple of other minor points which have nothing to do with the main idea. But 
the deep mystery is what I have described, and no one can go any deeper today". As 
Feynman mention in the lines above one could think that the interference effects are 
there because of the complex nature of the wave functions. Then it is natural to check 
what happens in the classical KvN case where, as we showed in the previous section, 
the wave functions have to be complex. We will actually show that, despite the com- 
plex nature of these wave functions, interference effects do not appear. This confirms, 
as Feynman suspected, that the mystery of quantum mechanics is deeper than that. 

If we want to describe a classical two-slit experiment we have to deal with a two 
dimensional problem. Let us call y the axis along which our beam propagates and x 
the orthogonal axis. We suppose that y = is the starting coordinate of our beam. 
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The centers of the two shts Ai and A2 are placed respectively at and — on a first 
plate which has coordinate t/p along the y axis. The final screen is placed at yg like in 
the figure below. 



XA 



i 



t Ai 



-XA 



A2 
Vf 



vs 



To simplify the problem we make the assumption^ that the motion of the particle 
along the y direction is known precisely. This means that at the initial time we know 
with absolute precision the position and the momentum of the particle, for example 



?/(0) = 0, p,(0) 



With this prescription we are sure that the beam will arrive at 



the two slits after a time tp = ypm/py and at the final screen after a time tg = ygm/py. 
In this way we can concentrate ourselves only on the behaviour of the particles along 
the X-axis. Suppose we consider, along x, a double-gaussian wave function with an 
arbitrary phase factor G{x,Px)'- 



ilj{x,px,t = 0) 



1 



V nab 



exp 



X 



2a? 



^ + ^G{x,p, 



(5.1) 



We assume a and h sufficiently large, i.e. the initial classical wave function sufficiently 
spread, in order to allow the beam to arrive at both slits. The evolution of the wave 



function will be via the free kernel of propagation ( |3.14| ) up to the time tp = yp-m/p^ 
that is when the beam arrives at the first plate. The wave function at the time tp will 
be: 



ip{x,p^,tf 



1 



exp 



1 



2a2 



X — 



VxVf 

Py 



A + ^g(x-^ V 



(5.2) 



Let us suppose that the width of the two slits is 26, then the particles which at time 
tp are outside of the two intervals Ai = {xa — S,Xa + S) and A2 = {—Xa — S, —Xa + S) 
are absorbed by the first plate and they don't arrive at the final screen at all. Using 
Feynman's words again: "All particles which miss the slit[sj are captured and removed 
from the experiment [l^"- Therefore the wave function just after tp can be rewritten 



in a compact way using a series of ^-Heavyside functions: 

i){x,p:,,tp + e) = N^p{x,pcc,tp) [Ci{x) + C2(x)] 



(5.3) 



similar assumption was made by Fcynman in his analysis of the diffraction experiment pj 
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where Ci{x) = 9{x — Xa + S)—9{x — Xa — S) is the function that parametrizes the sht Ai, 
C2{x) = 9{x + Xa + S) —9{x + Xa — S) is the one that parametrizes the sht A2 and is a 
suitable normalization factor chosen in such a way that: / dxdpx \ip{x,Px, tp + e)P = 1. 



Beyond the double slit we will propagate the ip of eq. (|5.3| ). With our choice of the 
cut-off functions Ci and C2, at + e the wave function ijj is different from only if 
X G Ai or ,T G A2. Since there is no limitation in the momentum along the x-axis 
we expect that the wave function ip will spread along x and, while time passes, it will 
become different from zero also outside the intervals Ai and A2. 

Using the kernel of evolution for free particles ( p. 141) we can obtain from (^.3|) the 
wave function at time tg = ysfn/p^ that is the time when the beam arrives at the final 
screen is: 



ipix.p, 



XI ts) 



N ■ exp ^ „ 
■{Ci{x - ap^) 



X 



PxVs 



262 



■exp 



iG 



X 



PxVs 



,Pa 



G2{x - apr,)} 



(5.4) 



where a = {i/g — yF)/Py- The probability density to find a particle in a certain point x 
on the last screen 



P(x) 



dpx\ip{x,px,ts)\' 



(5.5) 



We have to integrate over p^ because we are interested in the number of particles that 
arrive at the final plate independently of their momentum. At this point we notice a 
first important property: even starting from an initial wave function with an arbitrary 
phase factor G{x,Px), at time tg we have for the entire wave function a common phase 

factor of the form g(x — ^—y',Px\ see eq. (^.4|). So G will disappear completely 

^ Py 

in the evaluation of the modulus square and, consequently, in the P{x) of eq. (|5.5|). 
Therefore the phase G of the initial wave function ( |5.1|) cannot have any observable 
consequence in the figure on the final screen. 

The second important thing to notice is that, because of the properties of the 9- 
functions, we have that the cut-off term Gi + G2 in (|5.4|) is idempotent: 



(Ci + G2Y = Ci + C2 



(5.6) 



Therefore we can rewrite 



as: 



P(x) 



dp^\tp{x,px,ts)\'^ = N 



/oo 
dp^F'^{x,px,ts)Gi{x - dpx) + 
-00 



+ 



/oo 
dpxF'^{x,p^, ts)G2{x - dp^ 
-00 



where iV is a normalization factor and F is given by: 



F{x,p^,ts) = exp 



2a2 



X 



PxVs 

P'y 



pi_ 

262 



(5.7) 



(5.^ 
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Let us now re-arrange the arguments inside the ^-functions of the Ci and the C2 as 
follows: 



apx) 



C2{x — ap. 



0[-Px + 
o(-Px + 



X — Xa + S 

d 

X + Xa + S 
d 



0[-Px + 

o(-p.+ 



X + Xj, — 6 



(5.9) 



Remembering the properties of the 6'-Heavyside functions it is easy to realize that 

" X ^ S OC "-^A I S 

when P;j. is not in one of the two intervals: Di 

X + Xa — S X + Xa + S 



or D2 

a a 

there is no contribution to the modulus square. Therefore, 

a a J 

apart from the normalization coefficient N, we have that the final plot P{x) given by 
eq. (|5.7]) can be written as: 



Fix) 



dpx \'4){x,px,ts)\^ = N- 



dpx 



,3^5 Px-! ^s) + 



/ dpx 

JD2 



F'^{x,px,ts 



(5.10) 



where F is the function of eq. (|5.^ ). 

Now let us keep open only the first slit Ai and repeat the previous calculations. We 
can propagate the initial wave function ( |5.1|) up to the time tp when the system is 
again described by the ip{x,px,tp) of eq. ([5.2[). The difference is that now we have to 
parametrize only the first slit Ai. Therefore the second cut-off function C2 is identically 
zero. Since Ci itself is an idempotent function we can repeat the same steps as before, 
eqs. ( p.3| )-( pTT0D , freezing everywhere C2 to zero. As final result we obtain the following 
probability on the last screen: 



F{x) = k( dpxF\x,Px,ts 

JDi 



(5.11) 



where F is again given by eq. (|5.8| ). In the same manner keeping open only the slit 
A2 we will obtain that: 

P{x) = k( dpxF'{x,Px,ts) (5.12) 
ID2 



So, comparing ( |5.10| ) with ( |5.11| ) and ( |5.12| ), it is clear that when we keep open both 
slits Ai -|- A2 the total probability is the sum of the probabilities of having kept open 
first the slit Ai and then the slit A2. The first integral in ( |5.1CI| ) is then the probability 
for the particle to pass through the slit Ai, while the second integral is the probability 
to pass through the slit A2. So, even if we start from complex wave functions in the 
classical Hilbert space, every interference effect disappears. This is very clear from 
Figure |I| which shows the plot of the P{x) of eq. (|5.10|) with the particular numerical 
values Vs/Py = 2, a = b = 1,Xa = 1, S = 0.1. 

We will now perform a similar calculation at the quantum level and compare it 
with the previous classical experiment. In order to get an analytic result we will do 
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a simplification, that is we will assume that the motion along y is the same classical 
motion analyzed before. The reason for this assumption is that otherwise we would 
not be able to determine the time tp at which the wave function arrives on the plate 
with the two slits. Along x instead we will assume that the motion is fully quantum 
mechanical. So our overall approach to the quantum case is actually a " semiclassical" 
approach. Nevertheless this will be sufficient to see the difference with the purely 
classical case we have analyzed previously. 

In our semiclassical approach at the initial time the system along the ?/-axis is de- 
scribed by a double Dirac delta 6{y)6{py — Py) and this double Dirac delta evolves in 
time with the Liouvillian. In this way we know which is the time the particles arrive 
at the two slits. It is exactly the same as before. Along the other axis, x, we consider 
instead an initial wave function given by: 



ip{x) 



\/7ra 



exp 



X 



2a2 



(5.13) 



With this choice the mean value of both x and Px is zero at the initial time as in the 
classical case described by eq. ( |5.1| ). Making the above wave function evolve in time 
via the quantum Schrodinger operator, at time tp we obtain: 



i!{x,tt 



ma 



-exp 



mx 



2 ma^ + ihtf 



(5.14) 



Let us parametrize the two slits by means of the same series of 6'-Heavyside functions 
we have used in the classical case: 

Ci{x) = 9{x — Xa + S) — 9{x — Xa — S), C2{x) = 9{x + Xa + S) ~ 9{x + Xa ~ S) (5.15) 

Just after the wave function has passed the plate with the two slits we have that: 



ip{x,tp + e) = N ■ exp(-^- 



mx 



[Ci(x) + C2(x)] 



(5.16) 



2 ma? + ihtp, 

Using now the kernel of propagation ||l2i for a quantum free particle which is given by: 



K{X},,tb\Xa,ta) 



2mh{tb — ta 



m 



-1/2 



exp 



i m{Xb - Xg) 
2n{tb - Q 



(5.17) 



we get that at time tg the wave function is: 



ijix,ts) = m 



dxp exp 



im[x 



2(ma^ + ihtp) 



[Ci{Xp)+C2iXp)] 



where A^i is a new normalization constant. Differently from the classical case, the 
quantum kernel of propagation is not a simple Dirac delta and the previous integral 
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cannot be done explicitly. Anyway we can employ the properties of the ^-functions in 
order to rewrite (|5.18|) as: 



"XA+S 



im(x — XpY 



2h,{ts - tp) 2(ma2 + iht^ 

XA+5 r„-^^^ \2 .2 



+ / dxp exp 

J —XA—S 



im[x — Xp) mXp 



2h{ts - tp) 2(ma2 + iht. 



+ 

(5.19) 



In (|5.19|) we have two integrals of the same function over the two different intervals 
Ai = {xa — 5,Xa + 5) and A2 = {—Xa — 5, —Xa + 5). The results will be two complex 
numbers ipi and ip2 with different phases. So, differently from the classical case (|5.4|), 
the quantum wave function on the final screen ?/'(x, tg) has not a common phase factor 
and as a consequence the relative phases of ipi and il)2 will play a crucial role in giving 
interference effects. In fact if we re-write the final wave function as: 

V'(x, ts) = Ni IMx, ts) + Mx, ts)] (5.20) 

the probability on the last screen is given by the modulus square of tp{x,ts)'- 

P{x,ts) = 1^1 (a;,^^)^ + \Mx,ts)\^ + ri{x,ts)Mx,ts) + Mx,ts)r2{x,ts) (5.21) 

Note that the last two terms in the previous formula are not identically zero. If we 
make a plot of P{x, t^) as a function of x we see the evidence of the interference typical 
of quantum mechanics with a central maximum and a series of secondary maxima. 
This can be seen from Figure ^ which is the plot of P{x) in the case ts = 2,tp = 1, 
m = a = l,^ = l,5 = 0.1for two different distances of the slits: 2xa = 1 and 2xa = 2 
respectively. We can note the presence of six minima in the first case and of twelve 
minima in the second one. This in perfect agreement with the well-known relation 
that the distance Ax between two successive maxima or minima in an interference 
figure is inversely proportional to the distance 2xa between the slits. Therefore, even 
considering a quantum evolution only along the x-axis, the quantum wave functions 
create interference effects and the final result is a series of maxima and minima, like in 
the real experiment. 

Summarizing the results of this section we can say that if we make the evolution 
along the x axis with the Schrodinger Hamiltonian H, even starting from a real wave 
function, like ( |5.13| ), phases will appear during the evolution in a non trivial way and 



they will contribute to create interference effects. Instead in the evolution along x 
with the Liouvillian Ti, even starting from a complex wave function, like ( [5.1| ), the 
phase appears as a common factor for the entire ip on the final screen and so it does 
not contribute to {ipl"^ and it has not observable consequences. We feel that the two 
different behaviours are basically due to the different forms of the evolution operators 
in the classical and in the quantum case. 
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6 Conclusions 



In this paper we have shown, by means of simple and pedagogical examples, some of 
the differences between the operatorial approaches to classical and quantum mechanics. 
While in quantum mechanics the phases and the modulus of the wave functions are 
always coupled, eq. ( ^.331 ), in classical mechanics we can find a representation, the 
{q,p) one, where we can decouple them. Their being coupled in some representation 
of classical mechanics, like the (g, Ap) one, is only an apparent phenomenon. We think 
that this is the real profound feature of quantum mechanics which makes it different 
from classical mechanics: that in quantum mechanics there is no way to decouple the 
phases from the modulus by just going to a proper representation like we can do in 
classical mechanics . Of course we have not given a general proof of this in the sense 
that we have not really proven that in quantum mechanics this proper representation 



does not exist, but we strongly feel this is the case |13 



In this paper we have also performed the classical analog of the two-slit experiment 
and we have seen that interference effects do not appear. We have proved that by 
performing the detailed calculations giving the motion of the classical "wave functions" 
in the two-slit experiment. That was just a particular example and we do not know 
if it works in the same way in general like, for example, in the multiple-slit case or in 
other phenomena where at the quantum level there is interference. To do the detailed 
calculations in all these cases may turn out to be even more difficult than in the two- 
slit case. To by-pass those calculations we would like to find a general proof of the 
absence of interference effects in the KvN formalism. For sure in this general proof a 
crucial role will be played by the different form that the classical Ti, and the quantum 
H have. Another crucial role may be played by the universal symmetries present in 
classical mechanics and discovered in [Q. Those symmetries may be trigger a sort 
of superselection mechanism which may prevent the interference from appearing. All 
these features may be also those that, being absent in quantum mechanics, forbid the 
decoupling of the phase and the modulus. Work is in progress [O] on these problems. 



A Appendix 

In this appendix we will show that there is no contradiction in the postulate of KvN 
of having for ip the same evolution as for p = {ipl'^- The kernel P{ip,t\(pi,ti) for p is 
just ( p.9| ) which is a Dirac delta. So, postulating this to be the same as the kernel 
K{ip,t\ipi,ti) of propagation of the ip, in the case of a free particle we have: 

K{^,t\ipi,U = 0) = 6(q~qi- —)6{p-pi) (A.l) 
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Let us now use this expression to check what we get for the kernel of p knowing that 



= J dqidpi6(^q - Qi - ^^6{p - pi)ip*{qi,pi,0) ■ 

■ J dq'M5{q - q[ - ^^ip - p'Mq'.^p',, 0) (A.2) 
Now we can use the properties of the Dirac deltas to rewrite: 

p{t) = J dqidpidq[dp[5 {q-qt- 5{j) - Pi)S{pi - p'^) ■ 

■S(qi - q[ + {p^-p'i)-)r{q^.P^.Q)^{q'i.P^.Q) (A.3) 

V mj 



The integrals over the primed variables can be done explicitly 

t 



I dq-dp[5{pi-p'^5(^q^-q[ + {pi-p'i)-^iP*{qi,Pi,0)^P{q'i,p[^ (A.4) 
Substituting ( |A.4| ) into ( |A.3|) we have finally: 



pit) = J dqidpiK{ip,t\ipi,0)p{ipi,0) (A. 5) 

From this relation we get that the kernel of propagation of the p is the same as the 
one of the ijj and this proves that there is no contradiction in the KvN postulate. 
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Figure Caption 

Fig. 1: Classical Two-Slit Experiment. 
Fig. 2: Quantum Two-Slit Experiment. 
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